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An incommensurate spin density wave (Q phase) confined inside the superconducting state at high 
basal plane magnetic field is an unique property of the heavy fermion metal CeCoIns. The neutron 
scattering experiments and the theoretical studies point out that this state come out from the soft 
mode condensation of magnetic resonance excitations. We show that the fixation of direction of anti¬ 
ferromagnetic modulations by a magnetic field reported by Gerber et ah, Nat. Phys. 10 , 126 (2014) 
is explained by spin-orbit coupling. This result, obtained on the basis of quite general phenomeno¬ 
logical arguments, is supported by the microscopic derivation of the Xzz susceptibility dependence 
on the mutual orientation of the basal plane magnetic field and the direction of modulation of spin 
polarization in a multi-band metal. 

PACS numbers: 74.70.Tx, 75.30.Gw, 71.70.Ej, 74.25.Ha 


I. INTRODUCTION 


CeCoIns is a tetragonal, d-wave_^airing superconductor with the highest critical temperature Tc = 2.3K among 
all the heavy fermion compoundJiH^. The superconducting state of CeCoIns at a magnetic field above 9.8 T applied 
in the basal plane (Fig.l) of its tetragonal crystal structure co-exists with incommensurate antiferromagnetic (AF) 
ordering or spin density wave (SDW)p with = (0.45, ±0.45, 0.5) independent of the field magnitude. Its 2D 
incommensurate part q/c = (0.45, ±0.45,0) is parallel to the nodal directions of the d-wave order parameter 

A(k) = A(cos — cos ky). 

Here we use reciprocal lattice units. The existence of the magnetic order was first detected by the technique of 
NMRP and its precise field-dependence later determined^. The antiferromagnetic modulation is concentrated on 
the Cerium sites with amplitude m = O.WfiB {pb is the Bohr magneton) and polarized along the tetragonal axis. 
The incommensurate SDW is confined inside the superconducting phase, meaning that here superconductivity is an 
essential ingredient for SDW to develop. Different theoretical models have been proposed to explain why the SDW 
order occurs only in the high field-low temperature region of the d-wave superconducting stat^^®^. The choice between 
these models can be made with the help of results of neutron scattering. 

In the zero-field superconducting state, a spin resonance was found at a frequency Cli= 0.6 mev « 7K, which corre¬ 
sponds approximately to 3Tc. Initially, the resonance was observecP^ at a wave vector Q = (0.5,0.5, 0.5) associated 
with the nested parts of the Fermi surface corresponding to antiferromagnetic correlations. Theoretically Eremin and 
collaboratorJi^ have attributed the resonance to the proximity to the threshold of the particle-hole excitations con¬ 
tinuum which is at energy ui,. = min(|Ak| ± jAk+ql). Another scenario related to a magnon excitation was proposed 
by Chubukov and Gor’ko\ff^. 

Recent, more precise inelastic neutron scattering measurement^^ have demonstrated that the spin resonance is 
peaked at the same wave vector as the incommensurate static AF modulation in a high fields. Moreover, the fluc¬ 
tuations associated with the resonance are polarized along the tetragonal c-axis that corresponds to the direction 
of the ordered magnetic moments in the Q phase. Thus, the dynamical mode at zero field and the field induced 
static order share the same properties as if the resonance is a dynamical precursor of the Q phase. Also there were 
observed the Zeeman splitting of the resonance under magnetic field and the softening of the lowest energy mode of 
the Zeeman-split resonancd^lHaSI these observations point on the soft mode condensation caused by magnetic field 
as the mechanism for AF ordering formation. Such type theoretical model was put forward in Ref.16. 

In a Random Phase Approximation the spin susceptibility is 


x(q,w) 


1 - 


where w)) is the electron gas susceptibility in superconducting state and Uq is a momentum-dependent 

Hubbard-Coulomb repulsion potential. In this model the conditions for a collective excitation (called spin exci- 
ton) to occur are C/q3?ex^°)(q, w) = 1, and (q, w) ^ 1. The static antiferromagnetic state at some q is 

realized when the real part of the spin susceptibility along the tetragonal c-axis in the presence of a finite basal plane 
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magnetic field Hj^ exceeds the inverse constant of the AF interaction 

3?exi°;(u; = 0,q,i?j.)>C7-i. (1) 

In a two-dimensional model (Ref. 16) it was found that = 0, q/C; H±) in d-wave superconducting state increases 

with the field and exceeds the corresponding normal state susceptibility at fields essentially smaller than the param¬ 
agnetic limiting fielcP^. The physical reason for this behavior is that the incommensurate wave vector connects the 
points on the Fermi surface where A(k) = — A(k -|- q/c)- The same is true in real 3D case with modulation along 
Q/c- As a result, the tendency for AF instability in CeCoIns is much more effective in the superconducting d-wave 
state. 

An important observation made recentlj^^ is that the degeneracy between the two possible directions of antifer¬ 
romagnetic modulation Q/c = (0.45, ±0.45,0.5) is lifted by the magnetic field orientation. Namely, for field parallel 
to [110], a Bragg peak was with = (0.45,0.45,0.5) but no peaks corresponding to = (0.45,-0.45,0.5) were 
detected. For field precisely parallel to [100], the direction of the incommensurate part of AF modulation can take 
either of the two directions parallel to the gap nodes of dx^-y'^ pairing. But a tiny deviation of the field orientation 
from [100] toward [110] lifts this degeneracy and fixes the AF modulation along [the 110] direction. Correspondingly, 
a deviation of the field from [100] toward [110] fixes the AF modulation along [110] direction. In other words, the 
incommensurate AF modulation chooses that orientation which is the ’’most perpendicular” to the field direction (see 
Fig.2). 

The authors of Ref.25 have attributed this phenomenon to the presence in the Q phase of an additional modulated 
on atomic scale superconducting component with triplet p-pairing with order parameter Ag = i{AQcr)ay called pair 
density wave (PDW) interacting with d-wave order A and the AF order Mq 

Rcx*M§(AtdlQ-A(d^)t)+c.c. 

For the H j| [110] thep-wave state with maximal spin-susceptibility along field direction and the zeros of the order 
parameter along perpendicular to the field direction (dg oc {kx — ky)) does not disturb the d-wave superconducting 
state and can stimulate the emergence of Q-state with Q = (0.45,0.45,0.5). While possible phenomenologically this 
idea was not supported by any argument in favor of a specific mechanism for space modulated triplet pairing in this 
material. 

Another interpretation of the same phenomenon, developed quite recentlji^, is based on the AF modulation inter¬ 
action with the Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) modulation parallel to the magnetic field direction. One can 
remark, however, that the phase diagram for the coexistence of the superconducting Q-state and the FFLO state found 
theoretically in the previous paper of the same author^i^ does not resemble on the phase diagram of superconducting 
Q-state shown on Fig.l. Moreover, the isothermal measurement^^ at H jj [100] did not reveal an entropy increase 
at phase transition from the superconducting state mixed state to the superconducting Q-state which would indicate 
nodal quasiparticles in FFLO SC state. By contrast, a clear reduction of the entropy is found at a second-order 
transition at 10.4 T. This transition coincides with the incommensurate AF ordei®^. The observed reduction of DOS 
is in perfect agreement with the expectation for a SDW formation without the additional FFLO state. 

Here we propose an explanation of the Q-phase anisotropy governed by magnetic field based not on an imaginary 
additional ordering, but arising from the ordinary spin-orbit coupling. In the next Section, on the basis of quite 
general phenomenological arguments, we demonstrate, that in a tetragonal crystal under the basal plane magnetic 
field, the static magnetic susceptibility along the tetragonal axis xi? (q±,Hi) at finite q_L = {qx,qy) is largest either 
at qx II H_l or at q_L ± Hj_. As a result, if the maximum of the susceptibility occurs for a space modulation 
perpendicular to the field, then the inequality given by Eq.(l) is realized first on the Q/c direction, which is closer to 
being perpendicular to the magnetic-field direction. Thus, the degeneracy of directions of antiferromagnetic instability 
is lifted. 

The xi°i (w = 0, q, iLj^) calculated in a single band modeP^is completely independent of the mutual orientation of 
q and the basal plane magnetic field Hj^. One can show that this orientational independence persists in a single band 
metal even in the Abrikosov mixed state characterized by inhomogeneous superfluid velocity and field distributions. 
On the other hand, it is known that the spin-orbit interaction in a non-centrosymmetric tetragonal metal causes 
the magnetic susceptibility orthorhombic anisotropji^ Xxx — Xyy dx ~ dy- Such type anomalous susceptibility 
anisotropy in noncenrosymmetric CePtaSi has been mesured recently by polarized neutron scattering reported in 
Ref. 27. A similar phenomenon can be expected in a multi-band centrosymmetric material. CeCoIns is the multi-band 
metal that has been established by the de Haas-van Alphen measurement^^, by the band structure calculations and by 
the band spectroscopy studies^^"^. In the third section to illustrate the phenomenological conclusions we show that, 
owing to the interband spin-orbit interaction, the static spin susceptibility x^ (q_L,Hj^) in a tetragonal multi-band 
metal depends on the mutual orientation of the wave vector q_L and the magnetic field Hj^. The spin-orbit coupling 
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originating from the interaction of conducting electrons with the ionic crystal field can in principle be another source 
of violation of tetragonal symmetry by the basal plane magnetic field. 

The presented microscopic calculations of susceptibility are performed for a two-band tetragonal metal in the 
normal state. However, it should be stressed that antiferromagnetism and antiferromagnetic domain switching are 
phenomena originating from different mechanisms. As it was demonstrated in Ref. 16, the antiferromagnetism arises 
from an anomalous enhancement of the Xzz^(q/c) susceptibility in the d-wave superconducting state under the basal 
plane magnetic field. The domain switching is caused by spin-orbit coupling violating the crystal tetragonal symmetry 
at finite basal plane magnetic field and space modulation. The susceptibility Hj^) is proved to be dependent 

from the mutual orientation of the magnetic field and the direction of the space modulated spin polarization. This 
dependence takes place already in the CeCoIns normal state, but reveals itself in the d-wave superconducting state 
where the antiferromagnetic modulation developes. To demonstrate the violation of tetragonal symmetry by the basal 
plane magnetic field, it is sufficient to perform a microscopic derivation of the susceptibility in the normal state. The 
corresponding calculation in the superconducting state is much more cumbersome, but does not add any qualitative 
changes to the conclusions based on the normal state calculations. 


II. PHENOMENOLOGICAL APPROACH 


We consider a tetragonal paramagnet in a magnetic field having a constant basal plane (x, y) part and coordinate- 
dependent small additions 

H(r) = [H^ + SH^{r)]x + [Hy + SHy{r)]y + SH^{r)z. (2) 

Its free energy in the quadratic approximation has the following form 


T = J dV |a(Mj -F M^) + + 7 


\dx J \dy J 


fdmy dM^dMy 
dx dy 


+7-L2 


dM^ 

dx 




dM,i\ dM, 


dy 


dz 


+ <^1 dJ; 


dM, 


dM, 


dy dx 


+ S 2 {H, 


dM, 

dx 




dM, 

dy 


-HM 


( 3 ) 


Here, all terms besides the terms proportional to di and 82 have the tetragonal symmetry, whereas these two terms 
depend on the mutual orientation of the basal plane field and the direction of the space modulation of the M, 
component of magnetization. This dependence originates from the spin-orbit interaction. Disregarding the spin-orbit 
coupling corresponds to the equality (5i = 62 , recreating the functional tetragonal symmetry. For 5i> 82 ^ the direction 
of the M, basal-plane modulation is preferentially parallel to 


Hj. = HxX + Hyij. 

On the other hand, for 5i < 82 , the direction of the M, basal-plane modulation tends to be perpendicular to Hj^. 

One can formulate the same conclusion in terms of the susceptibility. Namely, by making the variation of Eq. @ 
with respect to the magnetization components, taking into account the expressions for the equilibrium parts of the 
magnetization 

Mx = iFx/So;, My = Hy/2a, 

and performing a Fourier transform, one arrives that the following equations for the magnetization response to the 
coordinate-dependent part of the magnetic field (5H(r): 


2 [a + ^{ql -f Qy) + ^,ql] SM^ + ^a:ykxky6My + ^x,qxqz5M, = SH^, (4) 
-fxykxkySMx; + 2 [u + J{q^ + qy) + J,ql] SMy + +-f±,qyq,5M, = SHy, ( 5 ) 
l±z{qxqzSMx + qyqzSMy) + 2[a,+ 'y{ql + q^) + 'y,ql + SiiH^qy - Hyqxf + 62 {Hxqx + Hyqy)"^] 6M, = SH,. ( 6 ) 
The solution of these equations yields the Fourier components of magnetization, 

6Mx{q) = XxxSHxiq) + XxySHy{q) + XxzSH,{q), 

6My{c[) = XyxSHxiq) + XyySHy{q) -f XyzSH,{q), 

6M,{q) = XzxSHxiq) + XzySHy{q) + XzzSH,{q). 


( 7 ) 

( 8 ) 
(9) 
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The coefficients 'yxy, and have relativistic smallness relative to the exchange-determined coefficients a and 7 . Ne¬ 
glecting the entanglement between the components of magnetization in Eqs. which gives only small corrections 

of order 0 {^ly) and 0 {'yxy'y±z) to the susceptibilities, we obtain 


Xzziq) = 


2 [ttz + '){ql + ql) + Izql + 5i {Hxqy - HyqxY + 62 [H^qx + HyqyY] 


Thus, for di > 82 the magnetic susceptibility along the tetragonal axis is largest for Hj_||qj_, where 


( 10 ) 


qj^ = {qx,qy)- 

On the other hand, for < 62 , the perpendicular mutual orientation of Hj^ and q_L corresponds to a maximum in 
the z component of susceptibility. This conclusion becomes evident if we rewrite the susceptibility as 

Xzz{q) 2 -h 7q5_-F 7 z9z + - <52)[Hj_ X q_L]2-H d2H^qi] ’ 

The derivation presented here demonstrates the dependence of the susceptibility on the mutual orientation of the 
field and the direction of modulation in the long wave limit. The effect, however, can be strong enough in the case of 
atomic scale antiferromagnetic orderings. 

Let us now show that the established Xzz susceptibility dependence from [Hj^ x qj^]^ really takes place in a two 
band tetragonal metal. 


III. MICROSCOPIC DERIVATION 


The Green function of a tetragonal two-band metal in an external magnetic field satisfies the equation 


HG=i 


with the 4x4 matrix Hamiltonian 


H = 


{iuj - ^i)ao -t- h_Lcr 

—Her (iuj 


Her 

^ 2 )cro + h_L cr 


( 12 ) 


(13) 


where 


^j(k) = £,(k) -/i, i = l,2 

are the band energies counted from the chemical potential, hj^ = fj,B is the Bohr magneton, and cr = 

(cto;, CTj,, CTj) are the Pauli matrices. The interband spin-orbit coupling^^ is given by the vector l(k) which is an even 

function 1(—k) = l(k) subordinating all the symmetry operations g of the tetragonal point group gl{g~^]i) = l(k). 

For the sake of concreteness, we can choose it to have the following form 

l(k) = l±kzikyX - kxV) + -i\\kxky{kl - kl)z, (14) 

where kx = kx/kp, ky = ky/kp such that 7 _l and qy have common dimensionality of inverse mass [1/m]. We are 
seeking a qj_ dependent z component of the spin susceptibility given by the equation 

Xi°i(q±) = -hIT 5] Tr ( ^ ) G(k -b q^/2,cc„) ( ^ ) G(k - q^/2,a;„). (15) 

Calculating the trace, we rewrite it in terms of the products of G matrix elements 


xi°)(q) = -/i|r^ E (-l)*+'G.,(k-bqj./2,cc„)G,,(k-qj_/2,a;„). 


(16) 
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In this complete expression, we seek the terms proportional to the combination [Hj_ x qj,]^. Let us take, for instance, 
the sum of the product of matrix elements 


[Gi 2 (k + q_L/ 2 ,a;n)G 2 i(k — qj_/2,a;„) + G 2 i(k + q_L/2,a;„)Gi 2 (k — qj_/2,a;„)]. 


(17) 


Among the many terms in this sum, we will keep only the Hj_ and q_L mutual orientation dependent terms. They are 


- Y, 




+ 


[{iuJn - 6)^ - - ll]k+q^/ 2 ^e[{h^ + ihyf{il^ + 

L>(k + qj_/2, a;„)D(k - q±/2, a;„) 

\{i(-Un — ^ 2 ) ~ ^_L ~ ^2]k-qj^/23?e[(/la; + iky) + ^i/)k+q^/2] 


where 

and 


D(k + q_L/2, w„)Zl(k - q±/2, a;„) 

77(k,a;„) = [(*a;„ - ^ - hi - + (h^l)i, 

Keeping the qj^ dependence only in 

we obtain 


^e[{hxihy) (*G + ^i/)k±qj^/2]’ 


(18) 

(19) 

( 20 ) 


U.,U)ri 


klljiuJn - 6(k))^ -h'j- lljk)] 

L> 2 (k,w„) 


{2[hj_ X q_L]2 - hiqi} . 


( 21 ) 


The sum over the Matsubara frequencies is easily calculated. However, the further integration over the Brillouin zone 
can only be performed numerically, taking into account the actual intraband quasiparticle spectra and the interband 
spin-orbital momentum dependencies. The knowledge of the quasiparticle density of states energy dependence near 
the Fermi surface is also crucially important. Analytically, we can write a rough estimation for orientation-dependent 
part of the susceptibility as 

2 


X. 


S“ oc ^>0 ( 






[H^ X q^]' 


( 22 ) 


Taking in mind that the modulus of the wave vector of antiferromagnetic modulation in CeColns is |q_L| « kp one 
can write the estimation of the absolute value of the susceptibility anisotropy as 

.2 AT I \ ( 7 -L^F^ ^23) 


anisl 

^zz 


2 AT (k-BH\ ('r±kl,\ 

(^j (^) 


CeCoIns is the heavy fermion compound with the electron effective mass about hundred times larger than the bare 
electron mas^^ m* « 100m. This means that the Fermi energy is quite small and can be of the order of spin-orbit 
interaction £p ~ 'y±kp. At the same time, the magnetic energy fipH in fields about 10 Tesla is about 10 Kelvin. 
Hence the anisotropy of susceptibility can have noticeable magnitude in comparison with the Pauli susceptibility 
Xp = 2No^%. 

Parametrically similar contributions to the orientation dependent part of the susceptibility originate from all the 
Green-function products in Eq. (16) with indices i ^ j. On the other hand the products with i = j give rise to 
contributions such as 


anis2 

Kzz 


(X h%Nq 


f fJ-BJWlpkjr 

{ el 


|H, X [>^1 (3!^ 


(24) 


Thus, we have demonstrated by direct microscopic calculation that the violation of tetragonal symmetry by a basal 
plane magnetic field discussed in previous section on a purely phenomenological basis, really takes place in multiband 
metals. 

Depending on the relative value and signs of expressions (22) and (241, the mutual orientations of qj_ and Hj^ can 
either be parallel or perpendicular to each other. The preferred mutual orientation can be changed at some pressure 
if the orientation-dependent part of the susceptibility changes sign. 
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IV. CONCLUSION 

Some time agcl^^, we have demonstrated the softening of spin resonance mode in the d-wave superconducting 
CeCoIns under a basal plane magnetic field at a wave vector q/c that connects the points of the Fermi surface with 
a finite gap A(k) = —A(k + q/c)- In the strong enough field this leads to the formation of static incommensurate 
AF state with two possible types of antiferromagnetic domains. Here, we have shown that the spin-orbit interaction 
in a tetragonal metal under the basal plane magnetic field acts to favor an inhomogeneous spin density modulation 
directed either perpendicular or parallel to the field direction. Hence, in general, only one type of antiferromagnetic 
domain is energetically favorable. This allows us to explain the puzzle of antiferromagnetic domain switching initiated 
by the basal plane magnetic field rotation observecP^in superconducting CeCoIns. 

In CeCoIns at ambient pressure, the incommensurate AF modulation prefers to be directed along the direction 
that is ”’he most perpendicular” to the field direction. At some pressure, the preferred orientation can change to ’’the 
most parrallel” one. 

The mechanism described for the orienting influence of the magnetic field on the direction of antiferromagnetic 
modulation has a general character and should reveal itself in an itinerant antiferromagnet. The phenomenon of 
antiferromagnetic domain switching is suppressed by domain pinning and can be observable only in clean enough 
metals, but not in doped antiferromagnets such as CeRhi_a;Co 2 ,In 5 . 

The susceptibility Xzz(<I-L) H-l) dependence from the mutual orientation q_L and Hj^ characterizing the intensity of 
spin-orbit coupling in a tetragonal material can be measured by neutron scattering. 

I hope that this paper will stimulate quantitative numerical calculations of Xzz(^±! H-l) both in the normal and in 
the superconducting states based on the real band structure of CeCoIns. 
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FIG. 1: Schemetical {H,T) phase diagram of CeCoIns. N and SC are the normal and the superconducting states corre¬ 
spondingly. The upper critical field in CeCoIns is mostly determined by paramagnetic limiting {Hc 2 {T = 0) ~ 11.7T) and due 
to this the phase transition to the superconducting state below T = QAT^ {T^ = 2.3 K) is of the first order (thick line on the 
figure(P. The Q-phase is the incommensurate antiferromagnetic state coexisting with the superconducting mixed state. 
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FIG. 2: (Color online). (Left) Basal plane magnetic field lies in the sector (0, ^) or in the sector (tt, ^), the AF modulation 
is directed along (1,1,0) direction. (Right) Basal plane magnetic field lies in the sector (f, tt) or in the sector (^, 27r), the AF 
modulation is directed along (1,1,0) direction. 






